function [c,ceq] = fminconstr1(x,eC,ell,theta,sigma,pi,rho,eta,i,M,A,B)

c = []; % no nonlinear inequality
ceq = sys1(x,eC,ell,theta,sigma,pi,rho,eta,i,M,A,B); % the fsolve objective is fmincon constraints